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Abstract 



' Given two simplicial complexes in TR d , and start and end vertices in each complex, 

we show how to compute curves (in each complex) between these vertices, such that 
the Frechet distance between these curves is minimized. As a polygonal curve is a 
complex, this generalizes the regular notion of weak Frechet distance between curves. 
We also generalize the algorithm to handle an input of k simplicial complexes. 
<y) . Using this new algorithm we can solve a slew of new problems, from computing 

a mean curve for a given collection of curves, to various motion planning problems. 
Additionally, we show that for the mean curve problem, when the k input curves are 
c-packed, one can (1 + e)-approximate the mean curve in near linear time, for fixed k 
and e. 

Additionally, we present an algorithm for computing the strong Frechet distance 
ly-j ■ between two curves, which is simpler than previous algorithms, and avoids using para- 

ph . metric search. 

o 

1 Introduction 

^ ■ The Frechet distance provides a way to measure the similarity between curves. Unlike the 

HausdorfT distance, which treats the curves as sets, the Frechet distance takes into account 
the structure of the curves, by requiring continuous reparameterizations of the curves. In- 
formally, the Frechet distance between two curves, it and a, is the minimum length leash 
needed to walk a dog when the person walks along it and the dog walks along a. 

In this paper, we are interested in extending this concept to facilitate solving more general 
motion planning problems. 
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Previous Work. The Frechet distance and its variants have been used to measure similar- 
ity between curves in applications such as dynamic time- warping [KP99], speech recognition 
[KH M + 98] . signature and handwriting recognition [MP99, SKB07J, matching of time series 
in databases |KKS05| . as well as geographic applications, such as map-matching of vehicle 
tracking data [BPSW051 IWSP06j . and moving objects analysis |BBG08al lBBG+08b| . 

Alt and Godau |AG95] showed how to compute the Frechet distance between two polyg- 
onal curves in IR d , of total complexity n, in 0(n 2 logn) tim^j]. It is an open problem to find 
a subquadratic algorithm for computing the Frechet distance for two curves. The decision 
problem (i.e., deciding whether the Frechet distance is smaller than a given value) has a lower 
bound of fi(nlogn) |BBK + 07j . Driemel et al. |DHW10] provided a (1 + ^-approximation 
for polygonal curves, that works in 0(N(e,ii,a) + N(l,n,a)\ogn) time, where N(e,n,a) 
is the relative free space complexity of two curves under simplification. In particular, their 
algorithm runs in 0(cn/e + cnlogn) time for c-packed curves. 

One can generalize the problem to consider an input of multiple curves. Dumitrescu 
and Rote |DR04] consider the problem of simultaneously minimizing the Frechet distance 
between all pairs of a set of k curves. They show one can get a 2- approximation in 0(n 2 log n) 
time, whereas the naive extension to k curves of the exact algorithm takes 0(n k logn) time. 
Buchin et al. |BBvK + 10] consider the problem of finding the median trajectory (polygonal 
curve) of a set of k trajectories in the plane that all share the same starting and ending points, 
where the median trajectory is a trajectory contained in the union of the k input trajectories 
that must cross at least half of the input trajectories in order to reach the unbounded face 
(though Frechet distance does not directly come in to play in their problem, it is related to 
the mean curve problem that we consider in Section 14. lj) . 

The notion of the Frechet distance can also be generalized to encompass distances between 
surfaces. Unfortunately, for general surfaces the decision problem is NP-hard [?]. In fact, 
whether the Frechet distance for general surfaces is computable is still an open problem. 
Recently Alt and Buchin [?] showed that the problem is semi-computable between surfaces, 
and polynomial time computable for the weak Frechet distance. The problem is hard even 
if the surfaces are well-behaved terrains, see Buchin et al. |BBS10j . 

Moving away from Frechet distances between surfaces, Alt et al. [AERW03] presented 
an 0(n 2 log 2 n) time algorithm to compute the Frechet distance between two graphs. Specif- 
ically, they require that one of the two graphs has to be entirely traversed and in the other 
graph we seek the path that minimizes the Frechet distance to the path of this traversal. 

Complexes. The notion of a complex (which is an abstract simplicial complex together 
with its realization), defined formally in Section |2T2| is a generalization of polygonal curves, 
triangulations, meshes, straight line graphs, etc. In particular, our algorithm uses complexes 
as inputs and as such would apply for all these different inputs in a verbatim fashion. 

Our Contribution. Given two complexes and start and end points in each one of them, 
we present a general algorithm that computes the two curves in these complexes that are 
closest to each other, under the weak Frechet distance, and connect the corresponding start 

1 In their paper, as well as ours, d is considered to be a constant. 
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and end points. In fact, the resulting curves that our algorithm produces are locally op- 
timal (see Remark 13. 5p . The running time of this new algorithm is 0(n 2 ) (where n is the 
complexity of the input complexes). Our algorithm can be interpreted as an extension of 
the algorithm of Alt and Godau [AG95J for computing the weak Frechet distance between 
polygonal curves. |AG95j and subsequent work on the Frechet distance has relied heavily on 
using the parametric space of the input complexes. Our main contribution is the usage of the 
product complex instead of the parametric space - this seemingly small change in viewpoint 
enables us to easily encode the, potentially very complicated, connectivity information of the 
input complexes in a simple way. This in turn allows us to observe that the classical Frechet 
distance techniques in fact apply to a wide array problems (many of which had not been 
previously considered) involving higher complexity input complexes. The notion of product 
complexes is implicit in previous work on this topic [AERW03| IBPSW05t IWSP06] . However, 
by bringing this concept to the forefront, we get a clean framework that facilitates our new 
applications. 

As concrete applications of our algorithm consider the following variants, all of them 
immediately solvable by our algorithm: 

(A) Frechet for paths with thickness. Imagine the classical setting of the Frechet dis- 
tance where a person walks a dog, but both the dog and the person might walk on paths 
that have non-zero width. That is, the input is two simple polygons (i.e., "thickened" 
paths) and one needs to compute the two paths of minimum Frechet distance between 
them that lie inside their respective polygons. 

(B) In a similar vane, consider a wiring problem: You are given a three dimensional model 
(of say a car or an airplane) specified by its mesh, and you are given a rough suggested 
path connecting two points in the mesh. Our algorithm can compute the optimal wiring 
path inside the model that is closest, under the Frechet distance, to the suggested rough 
path. 

Interestingly, this approach also extends to inputs of more than two complexes, and also 
to arbitrary convex functions between these different complexes. Specifically, consider a 
situation where the input includes k complexes Ci, . . . ,Ck- The reader might think about 
the complex Cj as the domain of the zth agent. Given a location in each of these complexes 
of their respective agent (i.e., a point Pi inside the complex Ci and the simplex A, C d that 
contains it) consider a scoring function /(pi, . . . , p&) that assigns a cost to the configuration 
(pi, . . . , pfc). Furthermore, assume that this scoring function is convex on the domain Ai x 
A 2 x • • • x A*;, and this holds for any combination of such simplices. Now, given that the 
agents want to move from some starting vertices Vi,...,Vk to ending vertices v[,...,v' k , 
the new algorithm can compute the synced motion of these k agents from the starting 
configuration to the ending configuration, such that the maximum cost of any configuration 
used throughout the motion is minimized. 

The reader might consider these settings a bit abstract, so here are a few examples of 
problems that can be solved using this framework: 

(PI) Mean curve. Given a set of k curves in IR d , find a new curve that minimizes the max- 
imum Frechet distance between this new curve and each of the input curves. Namely, 
this computes a mean curve for a given collection of curves. 
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(P2) One can compute the optimal way to walk k agents on k curves/complexes such that 
the maximum distance between any pair of agents, at any point in time, is minimized. 

(P3) Compute the optimal way for the k agents to walk on the k curves/complexes, such that 
the maximum average distance between any pair of agents is minimized (the average 
is over all pairs). 

(P4) Walk a pack of dogs while minimizing a weighted sum of the leash lengths (i.e. maybe 
some dogs need to be kept close since they like to chase squirrels). 

(P5) Motion minimizing the perimeter of the convex hull. Given k curves/complexes that k 
agents have to move on (in the plane), compute a motion from the start points to 
the end points, such that the maximum perimeter of the convex hull is minimized 
throughout the motion. 

The running time of all these algorithms for k input complexes of total complexity n 
is 0(n k ). In Section [5] we show that by making minor (realistic) assumptions about the 
input curves, for the median curve problem, one can remove the exponential dependence on 
k (however, the constant retains the exponential dependence on k). 

As a side problem, we also consider the problem when the input is two DAG complexes, 
which are directed acyclic straight line graphs embedded in IRA By considering the product 
space of two such complexes (instead of the parametric space) we show that the decision 
problem can be solved in 0(n 2 ) time. We then present a simple randomized technique to 
solve the general problem in 0(n 2 log n) time. In particular, this provides an alternative 
algorithm that computes the (strong) Frechet distance between two polygonal curves with- 
out using parametric search. Specifically, this algorithm is considerably simpler than the 
algorithm of Alt and Godau [AG95] . while matching its running time. Previous efforts to 
avoid the parametric search by using randomization resulted in algorithms that are slower 
by a logarithmic factor |vOV04t ?]. This new algorithm uses ideas applied for the problem 
of slope selection [Mat91j to the computation of the Frechet distance. See Theorem 16.31 for 
details. 

Organization. In Section [21 we define the Frechet distance and complexes formally, as 
well as introduce the key concept of using the product space instead of the parametric space, 
when defining the free space. Section [3] outlines the main algorithm of the paper, where 
it is shown that by applying the convexity property of the free space, our problem can be 
converted into the problem of computing the bottleneck shortest path. We also generalize the 
algorithm to handle k input complexes, as well as arbitrary convex functions. In Section H] 
we outline some applications of the main algorithm. In Section we show that when the 
k input curves are c-packed, one can solve the mean curve problem in near linear time. In 
Section [6], we present an algorithm for computing the monotone Frechet distance between 
two curves or between two DAG complexes. 
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2 Preliminaries 



2.1 Curves and the Frechet Distance 

Let Ti C M d be a curve; that is, a continuous mapping from [0, 1] to IR d . In the following, 
we will identify ti with its range 7r([0, 1]) C 1R d if it is clear from the context. 

A reparameterization is a continuous one-to-one function / : [0, 1] — > [0, 1], such that 
f(0) = and f(l) = 1. Given two reparameterizations / and g for two curves tt and a, 
respectively, define their width as 

width /iff (7r, a) = max \\rt(f(s)) - a(g(s))\\ . 

sG[0,l] 

This can be interpreted as the maximum length of a leash one needs to walk a dog, where 
the dog walks along tt according to /, while the handler walks along a according to g. In this 
analogy, the Frechet distance is the shortest possible leash admitting such a walk. Formally, 
given two curves tt and a in ]R d , the monotone Frechet distance between them is 

dj(n, a) = min widths Jtt, a) , 

/:[0,1H[0,1] h9K ' 

S :[0,l]->[0,1] 

where / and g are orientation-preserving reparameterizations of the curves tt and a, respec- 
tively. We will also be interested in the weak Frechet distance, where the reparameteriza- 
tions are required to be continuous but not necessarily bijections (i.e., one is allowed to walk 
backwards on their respective curve). In our problem we will be defining the curves in the 
respective domains. Hence finding curves that minimize the weak Frechet distance and find- 
ing curves that minimize the strong Frechet distance, are equivalent problems. The reader 
should note, however, that when the input domains are curves, our algorithm is equivalent 
to computing the weak Frechet distance between those curves. 

2.2 Complexes 

An n-dimensional simplex is the convex hull of n+1 affinely independent vertices. We call 
the convex hull of any m+1 vertex subset of the vertices of a simplex, an m-dimensional 
subcell (or face) of that simplex (note that a subcell is in fact an m-dimensional simplex). 
A proper subcell is one such that m < n. 

An abstract simplicial complex C\ = (P, 3 r ), is a set system. The elements of P are 
points and the elements of J are subsets of P called simplices . An abstract simplicial 
complex is downward closed; that is for any f 6 J, and T C \f, it holds that T 6 J. For 
our purposes, the ground set P will always be a subset of IR^. We also use the natural 
realization of the abstract simplicial complex (P, 3 r ), by mapping any simplex \& 6 H to 
rel(\I/) = CH(^), where CH(^) denotes the convex hull of Throughout our discussion we 
assume that for any \1/ e IF, we have = dim(C'H(\I')) + 1 (i.e. \l/ is affinely independent). 
We also require that our realization is locally consistent; that is V\l/, T 6 J, if ^ C\T $ 
then rel(^) fl rel(T) = rel($ n T). 
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Note, that the geometric realization of such an abstract simplicial complex does not 
induce a simplicial complex. For example, such an abstract simplicial complex might define 
a self intersecting polygonal curve, where two disjoint simplices ^ and T have that rel(\l/) 
and rel(T) intersect in their interior. In the following, we will refer to an abstract simplicial 
complex together with its realization as a complex. 

For a complex, C\, we will refer to any simplex in C\ as a cell of C\. The dimension of 
a complex is the maximum dimension of any of its cells. We say G C\ is a maximal cell 
of Ci if there is no Y G C\ such that ^ C T (note that a maximum cell is one such that 
dim(^) = dim{Ci)). 

A pair of simplices T are adjacent if $ C T or T C f . A simplicial path in a 
complex is a function : [0, 1] — > SF, such that: (A) For any A e J, we have that _1 (A) is 
a finite union of open intervals and points. (B) If </>(•) has only two distinct values (say A 
and on an interval [x, y] C [0, 1], then the simplices A and \1/ are adjacent. 

A curve tt C IR d parameterized over [0, 1] is a realization of a simplicial path 0, if for 
any t G [0, 1] we have that ir(t) G rel(0(t)) and cj){t) is the simplex of lowest dimension of 
3 that contains vr(t). In our applications, a maximal interval (x, y) such that (ft is constant 
corresponds to a straight segment of ir. In particular, when dealing with a curve tt C ]R d , 
we will assume that its associated simplicial path is also known. 

In the following we will abuse notation and refer to $ as a shorthand for rel(\I/). In 
particular, for a point p G IR^, we will say that p is in the simplex \I/ if p G rel^). 

2.3 Product Spaces 

Let C\ = (Pi, ^i) and C 2 = (P2, ^2) be two simplicial complexes in IR d . Consider the product 
space C\ x C 2 . Intuitively, we view the product space as a subset of the space JR 2d , where the 
first d coordinates are from C\ and the remaining d coordinates are from C 2 . With this view, 
Ci x C 2 is similar to a simplicial complex although the cells will be convex polyhedra instead 
of just simplices (in the literature this is known as a polyhedral complex). We define a 
cell (\l/, T) of Ci x C 2 to be the product of any cell \l/ from C\ with any cell T from C 2 . Its 
realization is the set rel(\I/, T) = rel(\l/) x rel(T). In the polyhedral complex C\ x C 2 , two 
cells (^, T) and (tf ', T') are adjacent if ^ is adjacent to ^' in d and T = T', or $ = ^' 
and T is adjacent to T' in C 2 . Also, note that C\ x C 2 is connected since, by assumption, the 
complexes C\ and C 2 are connected. 

Let 7r and a be curves with reparameterizations / and g, respectively. Let cell vr (-) and 
cell (T (-) be the simplicial paths associated with n(f) and u(g), respectively. Since the Carte- 
sian product of two continuous functions is continuous, we have that h(t) = (7r(/(i)), a(g(t))) 
defines a curve r = U t h(t) in C\ x C 2 , which we call the product curve of 7r(/) and cr(g). 
The curve r has a corresponding product cell path which is the function cell vr (7 (t) = 
(celL-^) , cello- (t)). (For the sake of simplicity of exposition, we are assuming here that 
cell vr (t) and ce\\ a (t) do not change their value simultaneously at the same time t.) 

For two complexes C\ and C 2 in IR d , and a parameter 5 > 0, consider a cell (Ax,A 2 ) 
in C\ x C 2 . For a point p = (pi,p 2 ) G (Ai,A 2 ), its elevation is the quantity elev(p) = 
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elev(pi, P2) = || Pi — P2 1| . The feasible region in the cell Ai x A 2 is the set 
F< j (A ll A 2 )={( I)! /)eI 2 ' i 



x £ rel(Ai) C R d 
y £ rel(A 2 ) C R d 
elev(pi, p 2 ) < 5 



The feasible region for C\ x C 2 (which we will refer to as the free spaced) is the set 

F< s {d,d) = U A F<s(A u A 2 ). 

Observation 2.1 Let tt and a be paths in C\ and C 2 , respectively, and let f and g be repa- 
rameterizations of n and a respectively, that realize the value 5 of the Frechet distance. 
The product curve, r, is contained in F<s(d,d)- Indeed, for any t £ [0,1], we have 
elev(jr(f(t)),a(g(t))) < 5, since f and g realize the Frechet distance between it and a. 

Observation 2.2 Consider a curve a in d x C 2 , such that for any point p £ a we have that 
elev(p) < 5. Then, the projection of this curve into the corresponding curves in d and C 2 
results in two curves Q\ and a 2 such that d^ci , a 2 ) < 5. 

Formally, for t £ [0, 1], let a(t) = (o"i(t), o~ 2 (t)) £ C\ x C 2 be a parameterization of a, and 
let cello- (t) = (celLj^) ,ce\\ a2 (t)) be its associated product cell path, such that for any t we 
have o~(t) £ rel(cell (J (t)). Clearly, o~i(t) and c 2 (t) are parameterized curves in the complexes 
C\ and C 2 , respectively. Furthermore, for any t £ [0,1], we have that Hcr^t) — a 2 (t) || = 
elev(cr(t)) < 5. As such, dg^Ui, a 2 ) < 5. 

2.4 Convexity of the Free Space in a Cell 

We need the following straightforward result. We include the proof for the sake of complete- 
ness. 

Lemma 2.3 Let F<s = F<s{Ci,C 2 ) be the free space of the complexes C\ and C 2 , both con- 
tained in lR d , for some 5 > 0. Then F<s(^, T) = -F<,5 n(rel(\I/) x rel(T)) is a convex set, for 
any cell T) of d x C 2 . 

Putting it differently, the elevation function elev(-) is convex over rel(\l/) x rel(T), for 
any cell T) of d x C 2 . 

Proof: Let ^ and T be simplices in d and C 2 , respectively, and let F = F<5(\1/,T). By 
the definition of free space, we know that F is just the sublevel set (i.e. the level set and 
everything less than that level) of the function h : 1R d x H d — > 1R, where h{u,v) = \\u — v\\, 
when applied to $ x T. It is known that the sublevel set of a convex function with a convex 
domain, is convex. Hence all we need to show is that h is convex (note that the domain is 
convex since \1/ and T are convex). 

So let u,u' £ ^ &ndv,v' £ T. We show that th(u,v)+(l-t)h(u'v') > h(t(u,v) + (1 -t)(u',v r )). 
for t £ [0, 1]. Equivalently, we show that the function h(t) = h(t(u,v) + (1 — t)(u r ,v')) is 
convex on the interval [0, 1], i.e. h(t) < (1 —t)h(0) + (t)h(l) (actually we need to prove such 



2 Note that here the free space is defined in terms of the product space, instead of the usual parametric 
space. 
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an inequality for all choices of u,u' G \I/ and v,v' G T, which holds since they were chosen 
arbitrarily) . 

Expanding out this function we get, 

h(t) = h(t(u,v) + (l-t)(u',v')) = h(v! + t(u - u'),v' + t(v - v')) 

= \\u' + t(u - u) - v ' - t(v - v') || = || {v! - v') + t(u + v ' - v! - V )|| . 

Hence h(t) is just the equation for the distance between a point on a linearly parameter- 
ized line and the origin. We have by Lemma 12.41 below that this function is convex and so 
we are done. ■ 

Lemma 2.4 The function representing the distance between a point on a linearly parame- 
terized line l(t) and the origin, is a convex function. Specifically, let a and b be vectors in 
]R d , then the function f(t) = \\a + tb\\, is convex. 

Proof: We know that f{t) is of the form, 



/(*) = , + thf = y/aP + pt + i, 



where a, (3, and 7 are some constants such that at 2 + (3t + 7 is non-negative. By the helper 
lemma below, however, we know such a function is convex. ■ 

Lemma 2.5 Consider the quadratic function at 2 + (3t + 7, where a, (3 and 7 are some 
constants such that the function is non-negative. Then, the function f(t) = \J at 2 + fit + 7 
is convex. 

Proof: Since at 2 + fit + 7 > for all t, it must be that a > 0, and the corresponding 
quadratic formula either has no roots, or a single root, which implies that j3 2 — 4«7 < 0. 
Now, 

= 2at + f3 = h(t)_ 
fU 2^at 2 + pt + 1 f(tV 
for h{t) = at + [3/2. Similarly, 

f{t)h\t) - f'(t)h(t) afjt) - {h{t)fl fU) = (f(t)) 2 - (Ht)) 2 /* 

1 1 ' (f(t)) 2 (f(t)) 2 (/(t))V« ' 

Now, since f(t) is always non-negative, we have that 

sign(/" (t)) = sign((/(t)) 2 - [hit)) 2 /a) = sign(at 2 + (3t + 7 - at 2 - f3t - (3 2 /Aa) 
= sign (7 - (3 2 /4a) = sign(4cry - f3 2 ) > 0, 

since a > and (3 2 — 4«7 < 0. ■ 
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2.5 Bottleneck Shortest Path Algorithm 



As a subroutine to our main algorithm, we need the following algorithm to compute bottle- 
neck shortest paths in linear time, which is accepted as folklore, but we include for sake of 
completeness. 

Let G = (V, E) be an undirected graph, with weight function w on the edges. For a given 
path p in G, let b{p) = max e€E ^w(e). We call a path connecting two vertices s and t in G 
an st path. We say that an st path, p', is a bottleneck shortest path, if among the set 
of all st paths, P, we have that b{p') = min pg p6(p). 

Lemma 2.6 Let G = (V, E) be an undirected graph, with weight function w on the edges. 
For any pair of vertices s,t £ V, one can compute in 0(\V\ + \E\) time the bottleneck shortest 
path between s and t. 

Proof: In the following, let p denote the bottleneck shortest path between s and t. First we 
check to see if s and t are connected in G and output oo if not. Otherwise, the algorithm 
proceeds in recursive stages. In each stage, we compute the median weight edge, e mec i, and 
consider the graph G< med = (V, E< med ), where E< med C E consists of all the edges of weight 
less than or equal to the median. We then compute the connected components of G< med . If 
s and t are in the same component of G< me( z then b(p) < w(e med ). In this case we recurse on 
the graph G< med since none of the larger weight edges can be in p. Otherwise, b{p) > w(e med ) 
and so we can contract each connected component of G< med down to a vertex, and recurse 
on the graph with these vertices and the edges in E >med — E \ E< med . In either case, at 
the end of a stage we remove all isolated vertices (this is also done before the first stage). 
Eventually, we will reach a stage with some small constant amount of edges in which case 
we can then solve the problem by brute force. 

Let 7T7.j and rij be the number of edges and vertices, respectively, in the ith stage. Com- 
puting the connected components, finding the median, and removing bad edges and isolated 
vertices all take 0(rij + J7ij) time (or better). Since we always remove isolated vertices at the 
end of a stage, we know that rii = 0(rrii), and hence each stage takes 0{rrii + rij) = 0{m,i) 
time. In each stage we delete (roughly) half the edges and so there are O(lgm) stages, and 
in each stage 0(rrii) = 0(2~ t m) work is done, and so 0(m+n) work is done in total (the n is 
included since the first check for isolated vertices will cost 0(n) time, and it might be that 
m = o(n)). m 

3 Computing Optimal Frechet Paths in Complexes 

We are given as input two complexes, C\ and C 2) along with corresponding start and end 
vertices Si, t 1 and s 2 , t 2 . We wish to compute the paths n and o in C\ and C 2 , respectively, 
that minimize the Frechet distance over all paths that start and end at the respective start 
and end vertices. 
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3.1 Algorithm 

We construct a graph G = (V,E), called the cell graph of C\ x C 2 . Specifically, each cell 
(Ai, A 2 ) of C\ x C 2 corresponds to a vertex t>(Ai,A 2 ) G V, an d for every pair i>(Ai,a 2 )> u (a' a 2 ) £ 
V we create an edge iff (A 1; A 2 ) and (A' 1; A 2 ) are adjacent in C\ x C 2 . For A x G Ci and A 2 e 
C 2 , the elevation of their corresponding vertex i> = i>(Ax,a 2 ) G ^ is elev(t>) = d(Ai,A 2 ), 
where d(Ai,A 2 ) = rmn pere i(Ai),qerei(A 2 ) elev(p, q) = min p6rd (Ai), q erei(A a ) ||p - q|| is the dis- 
tance between these simplices. The point realizing this minimum is the realization of the 
vertex v, and is denoted by rel(t>). 

The cell graph is clearly connected since C\ x C 2 is connected. As such, for any pair of 
vertices u, v G V(G) there exists a uv path in G. The elevation of a path p, denoted by 
elev(p), is the maximum elevation of any vertex in p. The lowest uv path in G is the uv 
path with minimum elevation. 

We compute the lowest st path in G (where s = (si,s 2 ) and t = (ti,t 2 )), in order 
to determine the desired curves with minimum Frechet distance. To this end, we set the 
elevation of any edge uv G E(G) to be elev(-uf) = max(elev(u) , elev(t>)). Clearly computing 
the lowest st path in a weighted graph is the same as computing the bottleneck shortest 
path, and so using the algorithm of Lemma 12.61 we can efficiently compute the lowest st 
path, p = v 1 . . . v m , in G. We return the polygonal path rel(t> 1 ) re\(v 2 ) ■ ■ ■ re\(v m ) C C\ x C 2 
as the desired curve (which by Observation 12.21 encodes the two desired curves and their 
reparameterizations) . 

3.2 Analysis 
3.2.1 Correctness 

As the following lemmas show, the cell graph captures the relevant information for our 
problem. 

Lemma 3.1 Let C\ and C 2 be two complexes, and let s\ and t\ be vertices of C\ and let s 2 
and ti be vertices of ' C%. Then, if there exists an s\t\ path it, in C\, and an s 2 t 2 path a, 
inCz, such that dj(7r, a) = 5 then there exists a V( sljS2 )V( tl ,t 2 ) path, p, in G(Ci,C 2 ) such that 
elev(p) < 5. 

Proof: Let / and g be the reparameterizations of 7r and a, respectively, that achieve the value 
S for the Frechet distance. By Observation 12.11 the product curve r = U t (Tt(f(t)), a(g(t))), 
defines a path in C\ x C 2 from (s%, s 2 ) to (ti, t 2 ) that is contained in the free space F<$(Ci, C 2 ). 
Let cell 7I - iCT (t) be the product cell path in C\ x C 2 that corresponds to r(t). Naturally, the 
value of cell 7I - )Cr (t) corresponds to a vertex in G, and let v(t) denote this vertex. It is easy to 
verify that the sequence of different vertices visited by v (t), as t increases from to 1, is a 
valid path in G. Indeed, a product cell path defines a sequence of adjacent cells of C\ x C 2 
as t increases from to 1, which corresponds to a path p — v 1 , . . . , v m in G. 
Observe, that for any t G [0, 1], we have that 

elev(v(t)) = eleven ^ = min ||p - q|| < ||7r(/(t)) - cr(g(t))\\ < 5. 

pGcelW(t), 
q£cella-(t) 
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As such, elev(p) = max, elev(V) = max t elev(v(t)) < 5. 



Lemma 3.2 Let C\ and C2 be two complexes, and let s\ and t\ be vertices of C\ and let S2 
and t2 be vertices of €2- Then, if there exists a v (si^Wti^) P a th P ^ n G(Ci,C2) such that 
elev(p) = 5 then there exists an sit\ path, it, in C\ and an Sit<i path, a , in C%, such that 
dy(7r,o-) = 5. 

Proof: Let p = v 1 . . . v m , where v 1 = V(si,s 2 ) an d v m = V(t lt t 2 )- Each vertex v % in p corresponds 
to a pair of cells A* = (A^, A 2 ), where A^ G C\ and A 2 G C2. Furthermore, for every z, there 
exists two points p* x G A l x and p 2 G A 2 , such that elev(p*) = \\p\ — p 2 || = d(A* 1 , A 2 ), where 

p*=(pi,pl)- " ^ V 

Observe, that for all the vertices of the polygonal path Z = p x p 2 . . . p' m , we have that 
elev(p l ) = d(A J 1 , A 2 ) = elev(V) < elev(p) = 5. As such, all the vertices of Z are in the free 
space F<s- 

For any i, the ith segment of Z is p*p* +1 . It corresponds to the edge v l v l+1 in the graph 
G, which connects adjacent cells in C± x C 2 - in particular, it must be that either A* C A l+1 
or A l+1 C A 1 . Assume the latter happens (the other case is handled in a symmetric fashion). 
We have that p l p J+1 C A*. Furthermore, by the convexity of the free space inside a single 
cell (i.e., Lemma [273]) . we have that p*p* +1 C A* D F<$. We conclude that Z C F<g. Since 
the two endpoints of Z are (s\, S2) = p 1 and (ti, t 2 ) = p m , Z corresponds to the desired paths 
Ti and a such that dg^vr, a) — 6. m 

Corollary 3.3 Let C\ and C2 be two complexes, and let S\ and t\ be vertices of C\ and let 
S2 and ti be vertices of €2- Moreover, let 7r and a be the paths in C\ and C2, respectively, 
that minimize the Frechet distance over all pairs of s\t\ and 52^2 paths. Then we have that 
d;j-(7r,cr) = 5 if and only if the lowest V( sliS2 )V( tltt2 ) path, p, in G(Ci,C 2 ) has elev(p) = 5. 

3.2.2 Running Time Analysis 

Computing the lowest st path takes + \E\) time by Lemma [2.61 Since a vertex in the 

cell graph represents a pair of simplices from C\ and C 2 , we know that |V(G)| = 0(|Ci||C 2 |). 
We also know that |£/(G)| = 0(|V^(G)|) since each cell in C\ x C 2 has at most 0(1) proper 
subcells (specifically 0(2 2d ) = 0(1)). Hence the running time of the algorithm is 0(n 2 ), 
where n = max(|Ci| , IC2I). 

Putting everything together, we get the following result. 

Theorem 3.4 Let C\ and C2 be two simplicial complexes, and n = max(|Ci|, IC2I). Given 
any pair of start and end vertices from C\ and any pair of start and end vertices from C 2 , we 
can compute, in 0(n 2 ) time, the paths n and a in C\ and C%, respectively, that minimize the 
Frechet distance over all paths that start and end at the respective start and end vertices. 

Remark 3.5 It is easy to verify that Theorem 13.41 yields a path that is locally as low 
as possible. Formally, if the solution in the polyhedral complex is a curve tt, then for any 
subcurve a C it, we have the property that for any other curve r, that has the same endpoints 
of a, it holds that elev(r) > elev(er). 
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When computing the Frechet distance for two curves for example, this property implies 
that the parameterization we get is never lazy - it always tries to be as tight as possible at 
any given point in time. 

3.2.3 Applications 

Frechet for paths with thickness. Given two polygons (maybe with holes) in the plane 
and start and end vertices in the two polygons, one can triangulate the two polygons and 
then feed them into Theorem 13.41 This results in two paths in the two triangulations that 
minimize the Frechet distance between the paths. As a concrete application, this can be 
used for solving the classical Frechet distance problem where the input curves have thickness 
associated with them and one can move in this enlarged region. Indeed, each "thickened" 
curve can be represented as a polygon, and hence we can apply the above algorithm. 

Wiring. The wiring problem, mentioned in the introduction, can be solved by immediate 
plug and play into the above result. 

Motion planning in planar environments. Consider the case where you need to plan 
the motion of two entities in a two dimensional environment, where they have to stay close 
together (i.e., Frechet distance) while complying with different constraints on which part of 
the environment they can travel on. As a concrete example, one entity might be a pedestrian 
and the other might be a vehicle. The pedestrian can not use the road, and the vehicle can 
not use the sidewalk or the parks available. Finding the best motion for the two entities is no 
more than solving the Frechet problem in this setting. Indeed, we compute a triangulation 
of the environment for the first entity, and then remove all triangles and edges that can 
not be used by the first entity. Similarly, we compute a triangulation for the second entity, 
removing the regions that are unusable for it. 

Now, applying the algorithm of Theorem l3.4l to these two triangulations (with the desired 
starting and ending points) results in the desired motion. 

Naturally, the algorithm of Theorem 13.41 can be applied in more general settings where 
the input is three dimensional, etc. 

3.3 Generalized Algorithm for k Complexes 

Let us recap the algorithm from the previous section. We considered finding the path in the 
product space (of two complexes) such that the maximum value of f(x, y) — \\x — y\\ among 
all the points (x, y) in the path is minimized. If we add an extra dimension for the value 
of /, then one can think of / as defining a terrain. Then the problem becomes computing 
a path that does not traverse high in this terrain. The free space was the sublevel set of 
/ for some parameter 5. Next, we defined the elevation of a vertex in the cell graph to be 
the minimum value of / for the cell that the vertex corresponds to. By observing that / 
was a convex function within each cell in the product space, we were able to argue that the 
value of the best path (i.e. lowest maximum value of /) was equivalent to the value of the 
bottleneck shortest path, and thus the problem was efficiently solvable. 
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With this abstract description, the only property of / that we used was that it was convex 
within each cell in the product space. Hence, we can conclude that the same procedure will 
work for any choice of /, so long as it is convex within each cell in the product space. 

We can generalize the problem even further. Earlier we considered only two complexes. 
However, there is no reason why we can not consider an input of k complexes, for some 
arbitrary integer k. In order to handle this case we generalize all our earlier definitions for 
two complexes in the following natural way. 

Let C\ = (Pi, 3^), . . ., Ck = (Pfe, Ujfc) be a set of k simplicial complexes in M. d . Consider the 
product space C\ x • • • x Intuitively, we view the product space as a subset of the space 
R kd . We define a cell (Ai, . . . , A fe ) of C = C\ x ■ ■ ■ x C k to be the product of k cells, where 
Aj e Ci, for % — 1, . . . , k. Its realization is the set rel(Ai, . . . , A&) = rel(Ai) x . . . x rel(Afe). 
In Ci x ... x Ck, two cells (Ai, . . . , A&) and (\I/i, . . . , ^k) are adjacent if there is a j such 
that for all % ^ j, Aj = and Aj is adjacent to ^ in Cj. 

We now are given a function / defined over lR kd that is convex for any cell rel(Ai, . . . , A&). 
As before, we build the cell graph G of the polyhedral complex C. Every vertex v of G 
corresponds to a cell A of C, and its elevation is the minimum value of / in this cell. 

As before, we are given start vertices Si, . . . , s k and end vertices t±, . . . ,t k in these k 
complexes. We compute the lowest elevation path between the vertex in G corresponding to 
(si, . . . , Sk) and the vertex in G corresponding to (t±, . . . , tk). Arguing as before, it is easy to 
show that the resulting path in the graph can be realized by a path in C that yields the k 
desired paths and their reparameterizations. As such, we get the following result. 

Theorem 3.6 We are given k simplicial complexes C\, . . . ,Ck, n = maxj \d\, start vertices 
Si G C\, . . . Sk G Ck, end vertices t\ G C\, . . . , tk £ Ck, and a function f : rel(C) — > 1R that is 
convex for any cell in the realization of C = C\ x ■ • • x C fc . 

Then, one can compute, inO(n k ^ time, k curves wi, . . . , 7r^ (and their reparameterizations 
■01, . . . , i[)k) connecting si,...,Sk to ti, . . . , t k , respectively, such that max ( f{^i{i^i{t)), • • • , 
7Tk{ipk(t))) is minimized, among all such curves and reparameterizations. 

4 Applications 

4.1 Mean Curve 

We are given k polygonal curves Hi, . . . ,7i k in IR d , and we would like to compute a curve a 
that minimizes the maximum weak Frechet distance between a and each one of the curves 

7Ti, . . . ,7Tjfc. 

For a set of points P C IR d , let r min (P) denote the radius of the minimum enclosing ball 
of P. 

Lemma 4.1 Let P(t) be a set of points in lR d moving linearly with t. Then, the function 
rmin(t) = r min (P(t)) is convex. 

Proof: Fix any three times, x < y < z, where y = ax + (1 — a)z for some a G (0, 1). Let 
Pi(t) denote the ith moving point of P(t). 
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Let y x (resp. v^) be the center of the minimum enclosing ball of P(x) (resp. P(z)), and 
let v(y) = av x + (1 — a)\i z . Observe that 

r min(y) = r min(P(«a ; + (1 — a ) z )) < max ||v(«X + (1 — a)z) — Pi(aX + (1 — a)z)\\ 

i 

< max(a ||v x - Pi(x)\\ + (1 - a) ||v 2 - Pi(z)\\) 

i 

= ar min (x) + (1 - a)r min (» , 

since the distance between a pair of linearly moving points is convex (for example by 
Lemma I2.4[) . ■ 

Using the lemma above, we get the following desired result. 

Lemma 4.2 Given k curves 7Ti, . . . , 7Tfc in IR d with total complexity n, one can compute, in 
0(n fc ) time, a curve a that minimizes maxj dj(7Tj, a), where d^(ni,a) is the weak Frechet 
distance between 7Tj and a. 

Proof: A cell in the polyhedral complex of 7Ti x • • • x -n k is the product of k segments (or 
points) in H d . For a point p = (p 1; . . . , p fc ) e 1R dk inside such a cell, consider the elevation 
of p to be /(p) = r m j n ({pi, . . . , p fc }). Lemma |4~T1 implies that /(•) is convex inside each such 
cell. As such, applying Theorem 13.61 to the given curves, using the function /(•), results 
in a parameterization that minimizes the maximum radius of the minimum enclosing ball 
throughout the motion. Since the center of the minimum enclosing ball (for continuously 
moving points) changes continuously over time, the curve formed by this center throughout 
the motion is a natural mean curve. Let a denote this curve. It is easy to prove that the 
maximum Frechet distance of a to any of the curves tti, . . . , iik is the minimum such value 
among all possible curves. ■ 

4.2 Walking a Pack of Dogs 

So suppose you have a pitbull, a chiwawa, a corgi, and a terrier. You want to walk all the 
dogs at the same time instead of walking each one individuallyll However, as before, long 
leashes are expensive, so you want to minimize the maximum length leash (among all the 
leashes) that you need to use. 

Formally, you are given k complexes, Ci, . . . ,Ck, and start and end vertices Sj,tj G Ci, 
for i = 1, . . . , k. The first complex corresponds to the person leading the dogs, and the 
complexes C 2 > • • • , Ck corresponds to the k — 1 given dogs. You wish to find the set of paths, 
7Tx, . . . , 7Tjfe, and corresponding reparameterizations, ^i, tp 2 , ■ ■ ■ , V'fcj sucn that, 



te[o,i] i>i 



max max ||7ri(-0i(t)) - Tti(ipi(t))\\ 



te[o,i] 



max /(7ri(V>i(t)), • • • , ir k (if> k (t))), 



is minimized, where f(p 



. . . , p fc ) = maxi || pj - pi 



3 Since clearly you are a person that is very concerned with efficiency. 
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Lemma 4.3 Given k polygonal curves 7ri,...,7Tfc of total complexity n, one can compute 
non-monotone reparameterizations of these curves such thatma.x t m&Xi \\iri(if)i(t)) — 7Ti(ipi(t))\\ 
is minimized. The running time of the algorithm is O (n fc ) . 

This works verbatim for complexes, and in this case the algorithm also computes the paths 
inside the complexes realizing the Frechet distance. 

Proof: We need to prove that the function /(pi, . . . , Pk) = maxj ||px — pj is convex within 
each cell in order to apply Theorem 13.61 

So, consider a cell A = (Ax, . . . , A&) G C = C\ x ■ ■ -xC fe . Its realized cell rel(A) = rel(Ai) x 
■ ■ • x rel(Afc) is a convex set. In particular, consider the functions of the form /»(pi, Pi) = 
|| Pi — Pi||, defined over rel(Ai) x rel(Aj), for 2 < i < k. Each of these functions are convex by 
Lemma [2731 on the domain rel(Ai) x rel(Aj). In particular, setting (?j(pi, . . . , Pfc) = /i(pi, Pi), 
for i = 1, . . . , k, results in k convex functions over rel(A). 

Clearly, /(pi, . . . , Pk) = maxj <7i(pi, . . . , p&), which is convex as the maximum of a set of 
convex functions is a convex function. As such, plugging this into Theorem 13.61 implies the 
result. ■ 

4.3 More General Settings 

From the previous example, consider the person and the dogs at any given time as vertices 
in space. The leashes are thus edges connecting the vertices. Hence in the above example 
the topology of the graph is that of star graphs (i.e. the person is at the center and the dogs 
are the ends of the star). The "weight" of each edge in the graph is the value of a convex 
function between the respective pair of vertices at a given instance of time (i.e. the distance 
of the person to a specific dog at a specific time). The general function we were trying to 
minimize was the maximum value over the functions between each pair of vertices. We were 
able to conclude that the overall function was convex because the maximum value of a set 
of convex functions, is a convex function. 

Let the above described graph be called a dependency graph. In general we can 
consider any topology for the dependency graph. More formally, between every pair of 
complexes we define a convex function (note that the zero function is convex, and so we can 
ignore certain pairs if we like). For our global function we can then take any function of 
these functions, which preserves convexity. For example, taking the maximum, the sum, or 
(positively) weighted sum of convex functions is again a convex function. Therefore, all of 
the applications (FH])-(F|4|) mentioned in the introduction are solvable immediately within 
this framework. 

4.3.1 Minimizing Perimeter of Motion 

We are given k complexes C\,...,Ck all with realizations in the plane. As before, we are 
given k starting vertices si, . . . ,s& and k ending vertices ti, . . . ,t&, in these k complexes, 
respectively. We are interested in computing the k polygonal paths (and their reparameter- 
izations) connecting these endpoints, such that the maximum perimeter is minimized. As 
before, to use the framework, we need to show that the perimeter function is convex inside 
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a cell of the resulting polyhedral complex. So, consider two points p = (pi, . . . , p^) and 
q = (q 1; . . . , q k ). We need to show that the perimeter function 

perim(t) = perim(tp + (1 - t)q) = perimeter (CH^tpi + (1 - t)qi, . . . ,tp k + (1 - t)q k })j 

is convex. This fact, which we state below as a lemma, is proved in [AClOj using the 
Cauchy-Crofton inequality. 

Lemma 4.4 ([AC10]) The perimeter of a set of linearly moving points in the plane is a 
convex function. 

This implies that the perimeter function is convex inside each cell of C — C± x • • • x C k , 
and hence the framework applies. We thus get the following result. 

Lemma 4.5 Given k complexes C±, . . .,C k all with realizations in the plane, and k starting 
vertices si, . . . , s k and ending vertices ti, . . . , t k , in these k complexes, respectively, then one 
can compute paths in these complexes, and their corresponding reparameterizations, such 
that the maximum perimeter of the moving points during this motion is minimized over all 
such motions. The running time of the algorithm is 0(n k ). 

The running time stated above is under the assumption that computing the minimum 
perimeter for k points whose locations are restricted by a cell of the polyhedral complex, can 
be done in constant time. This constant would depend on k, naturally. 

5 Computing the Mean Curve for c-packed Curves 

Driemel et al. [DHW10J introduced a realistic class of curves, called c-packed curves. We now 
show that when the k input curves are c-packed, one can compute a (1 + ^-approximation 
to the mean curve in O(nlogn) time, where 0() is used to emphasize that the constant 
depends on e and c, and exponentially on k and d (see Lemma 15.131 and Lemma 15.141 for 
details). This is a significant improvement over the algorithm for the general case, presented 
in Section H~Tl where the running time is 0(n h ). 

In this section, when we refer to the free space, it is meant with respect to the mean 
curve distance function. In particular, for k curves 7Ti, . . . , n k let d mean (7r 1 , . . . , n k ) denote 
the maximum distance of the mean curve to the tt^'s, for the optimum reparameterizations. 

5.1 Preliminaries 

5.1.1 Definitions and Lemmas 

We first cover the definitions and lemmas from [DHW10J that are relevant to our problem. 

Definition 5.1 For a parameter c > 0, a curve it in H d is c-packed if for any point q in 
H d and any radius r > 0, the total length of the portions of n inside the ball b(q,r) is at 
most cr. 
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Algorithm 5.2 Given a polygonal curve tt = qiq2q3 • • • and a parameter /i > 0, consider 
the following simplification algorithm: First mark the initial vertex q x and set it as the 
current vertex. Now scan the polygonal curve from the current vertex until it reaches the 
first vertex q, that is in distance at least \i from the current vertex. Mark qj and set it as 
the current vertex. Repeat this until reaching the final vertex of the curve, and also mark 
this final vertex. Consider the curve that connects only the marked vertices, in their order 
along 7r. We refer to the resulting curve tt' = simpl(7r, /x) as being the \i- simplification of 
7i. Note, that this simplification can be computed in linear time. 

We need the following useful facts about /^-simplifications from |DHW10] . 

Lemma 5.3 (i) For any curve tt in TR d , and fx > 0, we have that d^(7r, simpl(7r, //)) < /i. 
(ii) Let n be a c-packed curve in 1R d , let fi > be a parameter, and let it' = simpl(7r, //) be 
the simplified curve. Then, tt' is a Qc-packed curve. 

Observation 5.4 Let tt and a be two given curves, and let tt' and a' be their fi simplified 
curves, for some value fi. By Lemma \5.tA, dg?(7r, tt') < // and dj(o~,o~') < fi. Hence we have 
reparameterizations f and g such that ||7r(/(t)) — 7r'(t)|| < \i and \\a(g(t)) — o~'(t)\\ < \i for 
all t e [0,1] (without loss of generality we can assume these reparameterizations are bijec- 
tive). Let d^(ir',a') = 5, where d^(-, •) is the weak Frechet distance. Then we have that 
dj(7r,cr) < 5 + 2fi, since we can just map each pair (x,y) G (tt',o~') that is seen in the optimal 
(not necessarily injective) reparameterizations ofir' and a' , to the corresponding pair in (tt, a) 
determined by f and g. In particular, this implies that for curves n 1 , . . . , Tr k with correspond- 
ing ji simplifications 7r[, . . . ,ir' k , we have that d mean (^i, ■ ■ ■ ,^k) < d meon (7r 1 , . . . , ir' k ) + 2/i. 

Let TTi, . . . , 7Tfc be k given curves. The complexity of the reachable free space for these 
curves, for a distance 5, denoted by J^<s(^i, ■ ■ ■ , 7Tfc), is the total number of cells in the 
polyhedral complex with non-empty intersection with F<s(tti, . . . , ttu) such that there exists 
a path with elevation < 8 from the start vertex to that cell. 

Definition 5.5 For k curves tti, . . . , tt^, let 

N(e, 7Ti, . . . , 7Tfc) = max N<,j(simpl(7ri, s8) , . . . , simpl^, e<5)) 

be the maximum complexity of the reachable free space for the simplified curves. We refer 
to N(£, tti, . . . , Tik) as the e-relative free space complexityoi tti, . . . , tt^. 

5.1.2 Subroutines 

We now list the relevant subroutines from |DHW10] , which carry over directly for our prob- 
lem. 

Using the same procedure as in [D HWIOj . one can build a decider, decider(<5, e, tti, . . . , Tik) 
that runs in 0(N(e, m, . . . , iTk)) time (the only difference being that in our case the BFS ig- 
nores monotonicity) . Specifically, we have the following. 
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Lemma 5.6 Let m, . . . , 7r k be k polygonal curves in 1R d with total complexity n, and let 1 > 
e > and 5 > be two parameters. Then, there is an algorithm decider (<5, e, 7Ti, . . . , ir k ) that, 
in 0(N(e, n 1 , . . . , 7i k )) time, returns one of the following outputs: (i) a (1 +e)- approximation 
to d mea „(vri, . . . , 7r k ), (ii) d mea „(vri, ...,ir k )<5, or (iii) d mean (vri, ...,ir k )> 5. 

Definition 5.7 Given a finite set Z C IR, we say an interval [a, /3] is atomic if it is a 
maximal interval on the real line that does not contain any point of Z in its interior. 

Algorithm 5.8 For a set of numbers Z, let searchEvents(Z, e, 7Ti, . . . , tt^) denote the al- 
gorithm that performs a binary search over the values of Z, to compute the atomic interval 
of Z that contains d mean (7Ti, . . . , ir k ). This procedure would use decider (Lemma 15.61) to 
perform the decisions during the search. 

Lemma 5.9 Given a set P of n points in lR d , let be the set of all pairwise distances 
of points in P. Then, one can compute in O(nlogn) time a set Z of 0(n) numbers, such 
that for any y £ ( 2 ) ; there exist numbers x,x' £ Z such that x < y < x' < 2x. Let 
approxDistances(P) denote this algorithm. 

The following subroutine, from [DHW10] , will allow us to efficiently check intervals with 
bounded spread for d mean (7ii, . . . , n k ). 

Lemma 5.10 Given k curves tti, . . . , 7r k in lR d of total complexity n, a parameter 1 > e > 
0, and an interval [at, (3], one can compute a (1 + e)- approximation to d meo „(7Ti, . . . , 7Tfc) if 
^meani^i, ■ ■ ■ , ^k) £ [a, 0], or report that d mean (7r 1 , . . . , n k ) ^ [a, f3}. The algorithm, denoted 

by searchlnterval ([a, (3], e, 7T%, . . . , 7r k ), takes O^U(e,7Ti, . . . , ir k ) log ^f^Jtl^l \ time. 

We will also need a new subroutine, called solver, but first we prove the following easy 
lemma about MST's. 

Lemma 5.11 Let G be a graph with non-negative weights on its edges. For any two vertices 
u,v £ V(G), for the unique path r between u and v in the MST, we have that elev(r) < 
elev(cr), where a is any uv path in G, and elev(r) is the maximum weight edge along the path 

T. 



Proof: For the sake of simplicity of exposition assume that all the weights on the edges of 
G are distinct. 

Consider a uv path a in G. If a is contained in the MST then we are done. Otherwise, 
let e be any edge of a that is not contained in the MST. Introducing the edge e into the MST 
creates a cycle, where all the other edges on the cycle are lighter than e (otherwise, e must 
be in the MST). Therefore, we can replace e in a by the portion of this cycle connecting 
its endpoints. This new path a 1 has one less edge outside the MST, and it holds that 
elev(cr') < elev(cr). Continuing in this fashion, we end up with a path r' in the MST between 
u and v such that elev(r') < elev(<r). Since the path in the MST between u and v is unique, 
the claim now follows. ■ 
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aprxMean(e, 7Ti, ... ,717,.) 

(A) P = y(vr 1 )U---UV(7r fe ) 

(B) Z <— approxDistances(P) (Lemma 15. 9p . 

(C) [a, (3) <- searchEvents(Z, e, 7Ti, . . . , 7r fc ) (Algorithm 15.81) . 

(D) Call searchlnterval([a, 8a], e, 7Ti, . . . ,7Tfc) (Lemma 15. 10H . 

(E) Call searchInterval([/3/2, /?],£, 7Ti, 7Tfc). 

(F) 5^solver([2a,/3/2],7Ti,...,7r fc ) (Lemma EE]). 

(G) Return the value returned by searchInterval([<5/2, 35/2], e, 7Ti, . . . , ir^). 

Figure 1: The basic approximation algorithm. 

Lemma 5.12 Let tti, . . . , TTf. be k polygonal curves in IR d with total complexity n, 1 > e > 
be a given parameter, 5* = d OTeon (7Ti, . . . , 7Tk), and N = N(e, ttl, . . . , vr^). Let [a, (3] be 
an atomic interval that contains 5*, and such that for any //, // G [a,/3] ; simply, fi) = 
simpl(7Tj, /i') for i = 1, . . . , k. Then one can compute in 0(N log N) tzme, a value 5 such that 
5* G [5 — 2a, 5 + 2a]. Let this algorithm be denoted by solver([a, 0], tti, . . . , 7Ck) 

Proof: Let /i = 0. We run the algorithm of Lemma 14.21 on 7t[ = simpl^i, fi), . . . , i\' k = 
simpl(7Tfc, fi), except with the following modifications. First, instead of using the bottleneck 
shortest path algorithm of Lemma 1231 we will use Prim's algorithm, starting from the vertex 
that corresponds to the starting points of the curves, where we stop when we reach the vertex 
that corresponds to the ending points of the curves. Also, instead of explicitly computing the 
cell graph, we only compute the relevant parts of the cell graph on the fly as they are needed 
for Prim's algorithm. Note that if 5 is the elevation of the shortest path in the MST from 
s to t, then Prim's is guaranteed to stay within N<5(7i^, . . . , 7r' k ) until reaching t. Moreover, 
by Lemma 15.111 this path will be the lowest st path. 

This modified version of the algorithm computes a curve a that minimizes maxj dj(7r^, a), 
in 0(N log N) time, since we are running Prim's on an effective graph of size N (and where 
E(G) = 0(V(G))). Observe that since the /i simplification is constant on the interval [a,/3], 
5 is the same value that would be returned had we set fi = 5*. Also, again since the fi 
simplification is constant on this interval, by Observation 15.41 and considering fi — a, we 
know that 5* G [5 - 2a, 5 + 2a]. ■ 

5.2 Algorithm 

Given k curves, tti, . . . , 7Tfe, Figure [1] shows the algorithm to efficiently compute a 1 + e 
approximation to d mean (7r 1 , . . . , n k ). Note that the algorithm depicted in Figure [1] performs 
numerous calls to decider, with an approximation parameter e > 0. If any of these calls 
discovers the approximate distance, then the algorithm immediately stops and returns the 
approximation. As such, at any point in the execution of the algorithm, the assumption is 
that all previous calls to decider returned a direction where the optimal distance must lie. 
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5.3 Correctness and Running Time 

5.3.1 Correctness 

In order to apply the algorithm of Lemma 15.121 we first need to find an atomic interval (or 
sub interval) , [at, that contains 5* = d mean (7Ti, . . . , tt^), such that none of the \i simplifica- 
tions of any of the k curves change for any choice of /i G [a, 0\. Note that by the way in 
which jj, simplified curves are constructed, Algorithm I5.2[ if we consider increasing the value 
of \i from to oo, the only events at which any of the \i simplifications of any of the curves 
change, are when \i is equal to one of the distances between a pair of vertices on one of the 
curves. Hence if Y denotes the set of all pairwise distances between vertices in P (step ([A]) in 
the algorithm) then in order to apply Lemma 15.121 we want the atomic interval with respect 
to Y that contains 5*. Since it is costly to compute Y explicitly, we instead compute an 
0(n) sized set Z (step (IB]) ), such that each value in Y is 2- approximated by some value in Z. 
Step Q performs a binary search over Z, using decider, in order to find an atomic interval 
[a, (3] containing 5*. Since each value in Y is 2-approximated by some value in Z, we know 
that the interval [2a, (3/2] is a subinterval of an atomic interval of Y. Hence by Lemma [5.101 
we know that steps ([D]) and (jE|) ensure that [2a, (3 /2] is a subinterval of an atomic interval 
of Y that contains 5* (and if not, these steps returned a (1 + e)- approximation for 5*). By 
Lemma [5. 121 we know that in step (JFj) . when we call solver on the interval [2a, (3/2] we get 
a value 5 such that 5* G [5 — 4a, 5 + 4a]. However, flD]) guaranteed that 5 > 8a, since we 
checked the interval [a, 8a]. This implies a < 5/8 and so 5* G [5 — 4a, 5 + 4a] implies that 
5* G [5 — 5/2, 5 + 5/2]. Hence we have an interval with bounded spread which contains 5* 
and so by Lemma [5. 101 ( lG|) efficiently computes a (1 + e)- approximation for 5*. 

5.3.2 Running Time 

Let n = \P\ and N = N(e, 7Ti, . . . , tt^). By Lemma 15.91 that the call to approxDistances 
in ()B]) takes O(nlogn) time. Since searchEvents just preforms a binary search over the 
0(n) values returned by approxDistances by using decider, which runs in O(N) time by 
Lemma I5T6] we know that ([C]) takes O(Nlogn) time. Since [a, 8a], [(3/2,(3], and [5/2,35/2] 
are all intervals with bounded spread, we have by Lemma f5.101 that steps ([D]), (|E]), and 
f l(?]) run in 0(N log(l/e)) time. Finally, by Lemma 15. 12} the call to solver in line (EJ takes 
0(N log N) time. We thus have the following. 

Lemma 5.13 Let 7Ti, . . . , 7T& be k given polygonal c-packed curves in H d of total complexity 
n, let e > be a parameter, and let N = N(e, m, . . . , 7Tfc). Then one can compute, in 
0(Nlog(n/e) +nlogn) time, reparameterizations of the curves that 1 + e approximate the 
value of d mean (iii, . . . ,iik) . In particular, one can 1 + e approximate the mean curve of 

7Ti, . . . ,7l k . 

5.3.3 Free Space Complexity 

Lemma 5.14 For k c-packed curves 7Ti, . . . , ir^ in H d of total complexity n, and < e < 1, 
we have that N = N(e, 7Ti, . . . , n k ) = 0((c/e) k ~ l n). 
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Proof: Let 5 > be a fixed parameter, \i = e5, and 7i[ = simpl(7Ti, /x), . . . , tt^. = simpl^, /x). 

The free space in the polyhedral complex is partitioned into connected components. We 
must bound the size of the component which contains the start vertex, that is the reachable 
free space, H. Observe that one can charge a maximal dimensional cell in the polyhedral 
complex to an adjacent lower dimensional cell, since maximal cells contribute to N only if 
one of their adjacent proper subcells contributes to N. A non-maximal cell corresponds to 
some vertex v on one of the curves, and either a vertex or an edge from each one of the k — 1 
other curves. Consider a ball, b, of radius r = 25 centered at v. We now wish to count the 
number of features from the other curves (i.e. edges or vertices) that intersect this ball. 



To this end, consider one of the other curves, ti[. Let be 
the set of all features of that intersect b. Consider a ball, b', of 
radius 2r around v. Since r > /x and the edges of a /x simplified 
curve are of length > /x (with the exception of the last edge), 
every edge feature in Xi must contribute at least length /x to the 
intersection of b' and 7r| (note that if the feature is a vertex, then 
it is adjacent to an edge which contributes at least length /x). By 
Lemma I5.3[ the total length of n[ inside this b' is at most Ylcr. 
Therefore, 





Similarly, for each of the other k — 1 simplified curves, there are also 0(c/e) features close 
enough to v, that can be involved in a cell that contributes to N. Such a cell in the polyhedral 
complex involves choosing the vertex v, and one of these 0(c/e) features from each of the 
other k — 1 curves, and hence there are |Xi| • |Jf 2 | . . . • |-Xfc-i| = 0((c/e) k ~ 1 ) such cells. Since 
there are n vertices in total we thus have that N = 0((c/e) k ~ 1 n). ■ 

5.4 The Result 

Theorem 5.15 Let 7Ti, . . . ,7Tjt be k given polygonal c-packed curves in lR d with total com- 
plexity n, let e > be a parameter, and let N = N(e, 7Ti, . . . , n k ) = 0((c/e) fc_1 rx). Then one 
can compute, in 0(N log N) time, reparameterizations of the curves that (1 + e)- approximate 
the value of d mean (jri, . . . , 7Tfc). In particular, one can 1 + e approximate the mean curve of 
7Tx, . . . , vr fc in 0(N log N) = 0((c/s) k ~ l n\og n) time. 

6 Computing Optimal Frechet Paths for DAG Com- 
plexes 

In this section, we present a simple algorithm for computing exactly the monotone Frechet 
distance between two polygonal curves. This algorithm has running time 0(n 2 logn), and 
uses randomization instead of parametric search. In fact, the algorithm is considerably more 
general and applies to a wider class of inputs. 
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DAG complexes. Consider a directed acyclic graph (DAG) with vertices in K d , where a 
directed edge p — > q is realized by the segment pq. We refer to such a graph as being a DAG 
complex. Given two DAG complexes C\ and C 2 , start vertices S\ G V{Ci) , s 2 G V(C 2 ), and 
end vertices t\ G ^(Ci),^ G V(C 2 ), the problem is finding two directed polygonal paths 
7r 1 ,7r 2 in C\ and C 2 , respectively, such that: 

(A) The path 7Tj uses only edges that appear in C{, and it traverses them in the direction 
compliant with the orientation of the edges in d, for % — 1, 2. 

(B) The curve 7Tj connects Sj to t{ in Cj, for i — 1, 2. 

(C) The monotone Frechet distance between 7i"i and 7r 2 is minimized among all such curves. 

Note that this problem includes the problem of computing the monotone Frechet distance 
between two polygonal curves (i.e., orient the edges of the curves in the natural way and 
consider them to be DAG complexes). 

6.1 The Decision Procedure 

The algorithm is a direct extension of the work of [AG95J. Their algorithm relied on the fact 
that there was a clear topological ordering on the cells of the free space, and hence reachability 
information could be propagated. In this case, there is also a topological ordering (since it 
is a DAG). Hence, in the product space of two DAG complexes there is an ordering of the 
cells according to the underlying ordering of the two DAGs, and this ordering is acyclic. 

So, let C\ and C 2 be the two given DAG complexes, 5 a specified radius, and Si,s 2 ,£i,t 2 
the given vertices. The problem is to decide if there are paths between the start and end 
vertices in the corresponding complexes of Frechet distance at most 5. 

Algorithm. Compute the topological orderings of the cells (i.e., vertices and edges) of C\ 
and C 2 . In the resulting ordering it holds that A -<j A' if A appears before A' in this 
ordering, for % = 1,2, where A, A' G C;. 

We compute the product complex C = C\ x C 2 , and compute the topological ordering of 
the cells of C. Formally for A = (Ai, A 2 ) , A' = (Ai, A' 2 ) G C we have that A ^ A' if and 
only if Ai A[ and A 2 ;< 2 A' 2 . Clearly, the ordering -< over the cells of C is acyclic, and 
can be computed in linear time in the size of the complex. 

Now, just as in |AG95] . we start at the start vertex in the product space (s 1; s 2 ), visit cells 
according to their topological order, and compute the free space and propagate reachability 
information on the fly when we reach a new cell. 

Since we are working in the product space instead of in the parametric space, the two 
dimensional cells are parallelograms instead of squares. 

The reachability information is being propagated in a manner similar to |AG95] , except 
we propagate between adjacent cells, instead of neighboring two dimensional cells. Note, that 
no pair of two dimensional cells are directly adjacent, as there must be a one dimensional cell 
separating them. As such, for each edge (i.e., one dimensional cell) of C we maintain the set 
of reachable points. Unlike in |AG95] , the reachability information along a bounding edge in 
the product space might not be a single interval, since potentially multiple cells propagate 
to that bounding edge. However, by Lemma 12. 3[ we only need to compute the first point 
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(according to the ordering along this edge) that is reachable on this edge (notice, that an 
edge is always a product of a vertex of one curve and a directed edge of the other curve, and 
as such it has a natural ordering). 

In particular, when the algorithm visits a cell A in this ordering, it fetches all the cells 
that are adjacent to it and appear before it in the ordering. For each adjacent cell, the 
reachability information computed is of constant size, and hence we can compute the reach- 
ability information for the new cell in constant time. Indeed, the handling depends on the 
dimension of A: 

(A) dim(A) = 0: (A is a vertex), the algorithm computes if it is reachable from any of its 
direct ancestors, and if so we mark it as reached. 

(B) dim (A) = 1: (A is an edge), the algorithm computes the first point on the edge reachable 
from its direct ancestors. 

(C) dim(A) = 2: (A is a parallelogram), the algorithm uses the reachability information on 
the two incoming edges and the incoming vertex to compute the reachability inside the 
parallelogram. (Clipping the region to F<$ inside this cell.) 

As the algorithm visits the cells in a topological order, the work in maintaining the 
reachability information, can be charged to a cell's predecessors. As such, overall, the running 
time of the algorithm is linear in the complex size. 

The size of a DAG complex is the number of edges in it (since we assumed implicitly that 
the input DAG complexes are connected). Let n be the number of edges in the larger of the 
two DAG complexes under consideration. There are potentially 0(n 2 ) cells in the product 
space C. As such, the running time of the decision procedure is 0(n 2 ). 

Lemma 6.1 Let C\ and C2 be two DAG complexes, n be the number of edges in the larger 
of the two, si, t\ G Ci, S2, t 2 6 C2 be start and end vertices, and 5 > be a parameter. Then, 
one can decide, in 0(n 2 ) time, if there exists two paths tt% and ir 2 in C\ and C2, respectively, 
such that (i) 7Tj connects with t i} fori = 1,2, and (ii) d™^,^) < 5, where d™(-, •) is 
the monotone Frechet distance. Furthermore, if such paths exist, the algorithm returns them 
together with their respective reparameterizations realizing this distance. 

6.2 Using the Decision Procedure 

In the following, let C\ and C2 be the two DAG complexes under consideration. We outline a 
randomized algorithm to compute the value of the Frechet distance between the two curves 
in C\ and C2, that start and end at their respective start and end vertices, that minimize the 
Frechet distance. 

The algorithm needs to search over the critical values when the decision procedure changes 
its behavior. These critical values are the same as in Alt and Godau [AG95J (vertex-vertex, 
vertex-edge and monotonicity events). Indeed, for any pair of paths in the DAG complexes, 
the critical values for these two paths are the same as in [AG95J. As such, since DAG 
complexes are the union of paths, the critical values are the same. 

In the following, let 5* denote the actual minimum value of the Frechet distance. Given 
a parameter <5, let decider(<5) be the decision procedure described above. Let extract(a, b) 
be a procedure that returns all critical values determined by C\ and C2 whose radius is in the 
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compFr(Ci, C 2 , s x , s 2 , t u t 2 ): 

R: random sample of /i = An 2 critical values 
Sort R 

Perform a binary search over R using decider 
X = [a, b] <— Atomic interval of R containing 5* 
S <— extract (a, b) 

1/ S : all critical values in [a, b] 

Sort S 

x 4— Smallest value in S for which decider accepts 
// Computed using a binary search 

Return x 



Figure 2: The algorithm for computing the Frechet distance between two DAG complexes. 

interval [a, b). Suppose, for the time being, that the following subroutines have the following 
running times: 

(A) decider^) runs in 0(n 2 ) time (Lemma 16. ip . 

(B) extract (a, b) runs in 0(n 2 log n + k log n) time, where k is the number of critical values 
with radius in the interval [a, b]. 

(C) One can uniformly sample a critical value from the set of all critical values in 0(1) time 
per sample. 

The new algorithm is depicted in Figure [2j 
6.2.1 Computing the Critical Values in an Interval 

To complete the description of the algorithm, we need to describe how to implement ex- 
tract (a, b). For the interval X = [a, b], we need to compute all the critical values with radius 
in X. We can explicitly compute all the radii of vertex- vertex and vertex-edge events in this 
interval and sort them in 0(n 2 logn) time, where n is the number of edges (since there are 
0(n 2 ) such events in total and each radius can be computed in 0(1) time). Indeed, for a 
vertex- vertex event, its radius is the distance between the two vertices that define it. Simi- 
larly, the radius of a vertex-edge event is the distance between a vertex and an edge. Both 
types of radii can be computed in constant time, given the two elements that define them. 

In order to compute the radii of monotonicity events in X, we apply a variant of the 
standard line sweeping algorithm (i.e KDS). Specifically, for two DAG complexes C\ and C 2 , 
consider finding all monotonicity events between an edge e of C\, and pairs of vertices from 
V = V(C 2 ). To this end, place a sphere of radius 5 at each point of V with radius 5 = a. 
We now increase the radius 5 till it reaches b. The algorithm maintains an ordered list L of 
the intersections of the spheres with the edge e. The events in this growing process are: 

(A) The first time a sphere intersects e (this will create two intersections, if the intersection 
happens internally on e, since after this point the sphere will intersect e in two places). 

(B) When the intersection point of a sphere with e grows past an endpoint of e. 
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(C) When two different spheres intersect at the same point on e. At this point, the algorithm 
exchanges the order of these two intersections along e. The value of 6 when such an 
event happens is the radius of a monotonicity event. 
At any point in time, the algorithm maintains a heap of future events. Whenever a new 
intersection point is introduced, or two intersections change their order along e, the algorithm 
computes the next time of an event involving these intersections with the intersections next 
to them along e. 

It is clear that between such events the ordering of the intersections of the spheres with e 
does not change. Similarly, for a monotonicity event to happen on e, there must be a point 
in time in which the corresponding spheres are neighbors along e. Hence, this algorithm will 
correctly find all the monotonicity events. 

It takes 0((n + k) logn) time to compute all the relevant monotonicity events involving 
e and V, where k is the number of such events. We must do this for all edges of C\ and 
hence it takes 0((n 2 + k') logn) time to compute all the monotonicity events between edges 
of C\ and vertices of C2, where k' = ^ fcj and fcj is the number of monotonicity events in 
the interval [a, b] involving the ith edge of C%. Therefore it takes 0((n 2 + k") logn) time to 
compute all the relevant monotonicity events between C\ and C2, where k" is the number of 
such events (i.e. both those involving edges of C\ and those involving edges of C2). 



6.2.2 Sampling Critical Values 

We can uniformly sample critical values in 0(1) time, as follows. A vertex-edge event 
is determined by sampling a vertex and an edge, a vertex-vertex event is determined by 
sampling a pair of vertices, and a monotonicity event is determined by sampling a pair of 
vertices and an edge. Since we can easily uniformly sample vertices and edges in 0(1) time, 
we can therefore do so for critical events. In general, the decision of which type of critical 
event to sample would have to be weighted by the respective number of such events. 



6.3 Analysis 

Let R be the random sample of critical values, of size 0(n 2 ). The interval [a, b] computed 
by compFr contains 5*. The call to extract (a, b) takes 0(n 2 logn + fclogn) time, where k 
is the number of monotonicity events. The following lemma shows that k = 0(n 2 ). 

Lemma 6.2 Let X = [a, b] be the interval computed by compFr, and let c be some positive 
constant. Then, 

1 



Pr 



number of critical events in [a, b] > 2cnlnn 



< 
n c 



Proof: There are 2(™)n < n 3 possible monotonicity events, 2n 2 possible vertex-edge events, 
and n 2 possible vertex-vertex events. As such, the total number of critical events is bounded 
by Z = n 3 + 2n 2 + n 3 < 2n 3 . 

Consider the position of 5* on the real line. Let C be the set of the radii of all these 
critical events, and let U~ (resp. U + ) be the set of M = cn Inn values of C that are closest 
to 5* that are smaller (resp. larger) than it, and let U = U~ U U + . 
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If the number of values in C smaller than 5* is at most M, then there could be at most 
M critical values smaller than 5* in [a, b\. The same holds if the C contains less than M 
values larger than S*. As such, in the following, assume that both quantities are larger than 
M. 

The probability that the random sample R of size \x = 4n 2 picked by the algorithm, does 
not contain a point of U~ , is at most 

/ \U-\Y ( clnn\ 4n2 , n , N 1 

This also bounds the probability that R does not contain a value of U + . As such, with high 
probability, [a, b] contains only events in the set U. Namely, [a, b] contains the radii of at 
most \U~\ + \U + \ < 2M monotonicity events, with probability > 1 — l/n c . m 

Combining all our results, we thus have the following theorem. 

Theorem 6.3 For two DAG complexes, C\ and C 2 , of total complexity n, with start and end 
vertices Si,ti G Ci,s 2 ,t 2 G Ci, the algorithm compPr(Ci, C2, s±, ti, S2, t 2 ) returns two curves 
7Ti and 7r 2 , such that -k\ (resp. n 2 ) connects si (reps. s 2 ) to t\ (resp. t 2 ) in C\ (resp. C 2 ). 
Furthermore, the monotone Frechet distance between tti and ir 2 , is the minimum among all 
such curves. The running time of the algorithm is 0(n 2 \ogn), with probability > 1 — l/n c . 

Remark 6.4 The above result implies that given two polygonal curves in H d one can com- 
pute the Frechet distance between them, in 0(n 2 logn) time (this running time bound holds 
with high probability), by a simple algorithm that does not use parametric search. 

7 Conclusions 

In this paper, we showed that the algorithm for computing the (weak) Frechet distance 
between two curves can be extended to more general settings. This results in a slew of 
problems that can be solved using the new framework. 

Monotonicity. Our main algorithm from Section [3] is an extension of the algorithm of Alt 
and Godau |AG95] for the weak Frechet distance. It is natural to ask if the new framework 
can handle monotonicity. In Section[6], we offered a very restricted extension of our framework 
to this case, in the process presenting a new simpler algorithm for computing the monotone 
Frechet distance between polygonal curves. 

For more general settings, if the underlying complex is not one dimensional then it is 
not clear what monotonicity means. Even if we restrict ourselves to the case of k input 
curves, for k > 2, it is not immediately clear how to handle monotonicity efficiently, and 
we leave this as an open problem for further research. Interestingly, there are cases where 
monotonicity actually makes the problem easier. 
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Running Time. The running time of the general algorithm is 0(n k ) when handling k 
input complexes and is probably practical only for very small values of k. In Section [5] 
we showed that one can get a (1 + ^-approximation for the mean curve problem for k c- 
packed curves in 0(n log n) time. It should be possible to extend this same procedure to 
approximate, in a similar running time, some of the other problems that are solved by the 
general framework, under similar assumptions on the input. 
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